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The hot and dense strongly interacting matter created in collisions of heavy nuclei at RHIC 
energies is modeled with relativistic hydrodynamics, and the spectra of real and virtual photons 
produced at mid-rapidity in these events are calculated. Several different sources are considered, 
and their relative importance is compared. Specifically, we include jet fragmentation, jet-plasma 
interactions, the emission of radiation from the thermal medium and from primordial hard collisions. 
Our calculations consistently take into account jet energy loss, as evaluated in the AMY formal- 
ism. We obtain results for the spectra, the nuclear modification factor (R\ A ), and the azimuthal 
anisotropy (yj) that agree with the photon measurements performed by the PHENIX collaboration 
at RHIC. 

I. INTRODUCTION 

As they interact only electromagnetically with the surrounding matter, real and virtual photons have the potential 
to probe the detailed dynamical history of high energy heavy ion collisions. Their mean free path inside the hot and 
dense medium being much larger that its typical size, the photons will in principle leave the interacting zone without 
rescattering, reflecting directly the properties of the medium at the time they have been produced. The photon is thus 
expected to be a good probe for the quark-gluon plasma (QGP), the search of which has driven many experiments 
over the last years. We concentrate in this paper on the conditions that prevail at the Relativistic Heavy Ion Collider 
(RHIC) Q, and see if the experimental results obtained there are amenable to a theoretical interpretation in terms 
of new physics. 

Since the experimental detection of photons involves the entire collision, the QGP contribution might be hidden, 
or simply its effect reduced, by the sum of all other sources. It it thus essential to have robust calculations for 
those contributions, which includes the photons produced during the overlap of the nuclei (prompt contribution), the 
hadron gas contribution, as well as the background coming from the decay of mesons (^,77), after the thermal freeze- 
out. Considering that the background can experimentally be substracted, in principle by reconstructing the former 
mesons, we will concentrate on direct photons produced at RHIC in this paper. Particularly, recent studies 0, H, 0| 
have highlighted the role played by jets in real and virtual photon production. In Ref. @, calculations have suggested 
that the direct interaction of jets with the QGP would generate an inverse anisotropy, which can be traduced in term 
of a negative coefficient V2- However, those results were obtained using a longitudinal expanding QGP. 

In this work, the effect of the transverse expansion on photon production is evaluated. We use a 2D+1 hydrody- 
namical model, which has been applied recently with success to reproduce the characteristics of particle production 
at RHIC, such as momentum spectra, radial and elliptic flow. The spatial eccentricity in the model is also compatible 
with that inferred from experimental HBT measurements |6j. In SecJTTl the various sources of photons are presented, 
and the way the transverse flow enter into their expressions is shown. The Sec. IIIII presents the definition of all 
experimental observables that will be calculated, while in Sec. IIV1 the results are presented and analysed. Finally, 
Sec. |V] contains a summary and the conclusion. 



* Current address: Defence Research and Development Canada, 2459 Pie-XI Nord, Val-Belair, QC, Canada G3J 1X5 



2 




FIG. 1: Physical processes without collinear effects. 

II. PHOTON PRODUCTION 
A. Jet-thermal processes 

From finite-temperature field theory Q , the production rate of virtual photons with momentum p, invariant mass 
M and energy E is 

F dW\ 1 Tr[lmn R ] 

where Tr [lmll fl ] = ImTJ^ M is the trace of the imaginary part of the retarded photon self-energy. In the hard thermal 
loop (HTL) resummation formalism , the non-collinear processes contribution up to next to leading order in g s are 
shown in Fig. [T] The filled circles in this figure indicate resummed propagators. From relativistic kinetic theory the 
production rate of those non-collinear processes induced by jets have the following form: 



d 3 p 



= J b ) E d^d t ■ ( 2 ) 



The color and spin degeneracy factors are respectively iV c = 3 and N s = 2. The phase-space distribion of incoming 
jets, initially created at position rj_, assuming a Bjorken 77 — y correlation [l(| is 



6tqt cPgrdy 



4 c + Vx, q , t, b) = ^"L? J A 2 s {r > - y) » (3) 



where r = \A 2 — z 2 is the proper time, the z-axis being the beam direction, and 77 = l/21n[(f + z)/(t — z)] is the 
space-time rapidity. The initial profile of jets in the transverse plane is obtained by 

Tab{o) 

where T4 and Tab are the thickness and overlap functions, which are evaluated with a realistic Woods-Saxon 
distribution [Tlj . The initial momentum distribution of jets, at a given impact parameter b, is obtained by 



dN iet (Q,b) 



cPqxdy 



= T AB (b) dx a g a {x a , 

t=0 a,Kc J 



Q)g b (x b , Q) 



da a+b^c+jet 2XaXb 

X dt ir{2x a -2^ev) () 

where y/s — 200 GeV. Isospin effects are included in the parton distribution function (pdf) by 

! Z A — Z \ 
9a{x a ,Q) = ( -rfa(x a , Q) H -j — f a * [x a , Q) J R(x a , Q) , (6) 

where f a is the parton distribution function inside proton [T3 |. The second term, corresponding to pdf inside the 
neutron, is obtained by the following substitution of parton species: a* = (d, u, d, u, s, s, g) for a — (u, d, u, c?, s, s, g). 
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Shadowing effects are inciuded in the function R{x„. : Q) [i3j |. The factorization scale Q is assumed to be qr- We use 
a NLO factor K = 1.7, which, according to Ref. is almost (^-independent. The jet distributions evolve in time 
according to [HI 
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where the k integrals run from — oo to oo. The transition rates in the laboratory frame are 



dV q qa {q,k) 
dkdt 
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^jet 



dT q qg (q ,k ) 



dkodto 



(8) 



where dT q qg / dkodto are evaluated in the fluid local frame moving with a velocity /3 relatively to the laboratory frame, 
and (1 — vjet • /3) represents the Jacobian of the dtodko — > dtdk transformation. The transition rates in the local 
thermal frame are taken from the AMY [3] formalism, and include radiative energy loss through gluon bremsstrahlung 
with LPM effect. Processes like qq annihilation and absorption of thermal gluons are also included in the model. 
The strength of the transition rates is controlled by the str ong coupling constant a s and the temperature T. The 
temperature dependence of a s is obtained from lattice QCD [16j |. As the jets propagate in the QGP, the parameters 
(3 and T which depend on the position of the jets and the time, are directly extracted from the hydro model. All 
jet-medium interactions cease when the critical temperature T c is reached. The jets distribution, appearing in Eq. ([2]), 
is evaluated in the laboratory frame, so must be the quark to photon transition rates. Since E /d 3 po is a Lorentz 
invariant, the Jacobian for the Lorentz transformation is simply dto/dt = \/l — |/3| 2 - We thus get 



E 



dr -n 7 -coll(S>P) = ^WEo 



<rn 7 -cou(90,P0,M) 



(9) 



d 3 pdt d 3 podto 

The quark to photon transition rates in the local thermal frame have been calculated in Ref. 0]. The photon's 
momentum and energy in that frame are 
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(10) 



The jet's momentum follows the same transformation rule, but with E q = q since we assume the jets to be massless 
(I Vjet = 1|)- The yield of virtual photons produced in non-collinear processes induced by jets in the expanding medium 



E dN ™?- mn (M) 



d 3 p 



d A xE 



dW 



(M) 



d 3 p 

drr drjd 2 x± N C N S 



d 3 



q 



(27T) 



3 b ) \/1-P(t,v,x±) 2 Eo 



al noil -coll 

d 3 p dt 



(11) 



The retarded photon self-energy 11^ M for collinear processes with LPM effects, i.e. including an infinite sum of 
diagrams with different number of scatterings with soft gluons (see Fig. [2]) , has been extended from real photons Q 
to virtual photons in Ref. [T3|, in the limit M <C E. The dilepton production induced by jets in collinear processes 
has finally been calculated in Ref. [T^ for a longitudinal expanding QGP. The quark to photon transition rates for 
collinear processes have been extracted from il^ using Eq. (TTJ) and relativistic kinetic theory: 



EQ~J3 — 



= 12 



dqo 



qpEp 
(2n) 3 p 2 



nFD(qo, 



^ r coTr(9o,Po,M) 



dpodto 



1 - 



0(po - qo) 



(12) 
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FIG. 2: Bremsstrahlung and annihilation processes with LPM effect. 



where npD is the Fermi-Dirac phase-space distribution function. The 0(po — qo) function is included to avoid double 
counting in the annihilation process, since ". fd ((Zo)"-fd(po — Qo) = ^fdCpo — < Zo)"-fd( ( ?o)- After substracting from 
^^coU 7 ^ ne l ea ding order annihilation (Fig. QJ,), again to avoid double counting, the photon yield from the collinear 
processes (f2j et ~ 7 and ?/j t ~ y 7 = y for M <gi px) induced by jets in an expanding medium is finally 



d 6 p J J J P 2 d 2 q T dy 



dp dt 



The Jacobian of the dto dpo — > dt dp transformation is 

Ei ( i 

Po \E p 



(13) 



J = ^(l- P -^\- (14) 



B. Thermally induced processes 



The photons produced during the thermal phase due to the collisions of thermal particles, are calculated simply by 



E dN ^t ima} (M) = 
d A p 



d 4 xE, 



dW 

ul "-thermal 

d 3 p 



(15) 



where the photon production rate is evaluated in the local thermal frame. The QGP induced processes also correspond 
to the diagrams shown in Figs. Q] and [H with the difference that the incoming particles are now thermal partons 
rather that jets. The production rates for the QGP induced photons are from Refs. [H an d fl7| , while the hadronic 
gas (HG) production rates (mesonic as well as baryonic) used here are introduced in Ref. [19|. Defining /qgp as the 
QGP content in the mixed phase, the thermally induced radiation is 
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where T^ ec denotes the decoupling temperature. 
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FIG. 3: (Color online) Distribution of real prompt photons in p — p collisions at RHIC. Data points are from PHENIX [25| . 



C. Prompt photon production 



The prompt photon production in p — p collisions is given by 
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(17) 



q T =p T /z 



Here, the renormalization scale has been implicitly set equal to the factorization scale Q and the leading order 
expression for the cross-section dcr a +b^ c + 7 * and da a+b ^ c+ d can be found in Refs [20( and [2l|. The iiT-factors include 
NLO effects. We assume, for M -C pr, that K(M,p,Q,Qp) K(p,Q,Qp), and we evalutate them using the 
numerical program from Aurenche ef aZ. |2^ |. The real photon vacuum fragmentation function comes from Ref. (23j 
while for M ^=0, we take the leading order result The splitting between the direct and fragmentation contributions 
is arbitrary and depends on the choice of fragmentation scale Qf ■ So, only the sum of the two contributions has a 
clear interpretation. We nevertheless set Q = Qf, and ajust Qf to fit the recent PHENIX prompt (real) photon 
results in p — p 

Using Q = p T j \/2, we obtain a nice agreement with data, as shown in Fig. [3l 
In nucleus- nucleus collisions, the direct contribution is simply given by 



E^±(M,Q,Q F ,b) =T AB (b) 



d 2 pxdy 



(M, Q, Q F ) 



(18) 



where the pdf functions f a are replaced by g a , defined in Eq. ([B]), to include isospin and shadowing effects. The frag- 
mentation contribution, like the jet-fragmentation contribution in pion production, could suffer high-px suppression 
due to the quenching of jets in the QGP. Remember that the distribution of jets has been calculated in Eq. §5§ by 
setting Q = qx, while in Eq. (jTTJ) , the scale has been set to Q = pt/V% ■ Thus, in order to make a connection with 
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the distribution of jets, we approximate the fragmentation contribution in nucleus-nucleus collisions as 
dNZl „ r . r 1 dM t (d,Q' = g T ,6) 



frae 



(M,Q,Q F ,b) ~ J2 J <Prj_V(r x ) J d lT 



Pt 



d 2 q T dy 



l=qq,g 

xH{q T , Q , Q) D r/t (M, z = p T /g T , Q F ) . 



(19) 



The effective function H(qx 1 Q , Q) has been introduced so that if the jet is not quenched in the medium, we get the 
following: 



E 



dNL 
d 3 p 



■(M,Q,Q F ,b) 



no E—loss 



d 2 p T dy 



(20) 



For Q = Q F = Pt/V? and Q' = q F , we obtain that H(qT,Q ,Q) ~ 1.9. The quenching of jets in the medium will 
depend on the path lenght d, which depend on the propagating direction of jets and the position r_L in the transverse 
plane where they have been created. 



III. DILEPTON AND REAL PHOTON PRODUCTION 



The total direct photon yield is the sum of all contributions discussed in the preceding sections: 

^£aL(M, b) = E ^" 11 + + thermal + E <^t + > ^ 

d 2 pTdy ' d 3 p d 3 p d?p d 3 p d 3 p 

where d 2 p F = pr dpr d<fi. The real photon spectrum is simply obtained by 

(b) = ^%1. (M = 0, b) , (22) 
a^prdy a^prdy 

while the dilepton spectrum is 



dM 2 d 2 p T dy^ * = ^d^rfy^ P ^ * Vcut, Pt, M) . (23) 

The multiplicative factor P(|y e ±| < y CVL t), defined in Ref. [H, is introduced in order to take care of geometrical 
acceptance of any detector. Finally, on top of the photon yield, the nuclear modification factor and the azimuthal 
anisotropy coefficient provide important informations about behaviour of jets in the medium. They are respectively 
defined by 

p7 v C d^dN^b)ld 2 p T dy 

2TrT AB {b)da^ ompt /d 2 p T dy 

and 

ft* #cos2^ dN^(b)/d 2 PT dy 

v4(b,PT,V) = — — o . (25) 

2 f^d0dNt(b)/d 2 pTdy 

In this study, we examine only the production of mid-rapidity photons, such that y = is set in all the above equations. 



IV. RESULTS 



The thermal fluid dynamical evolution is described by by the longitudinally boost-invariant (2+l)-dimcnsional 
hydrodynamic code AZHYDRO [111, H3 with the EOS Q equation of state which matches a non-interacting QGP 
above T c to a chemically equilibrated hadron resonance gas below T c at the critical temperature T c = 164 MeV. We 
assume an early start of the hydrodynamic evolution at Ti=0.2fm/c in order to be able to account schematically for 
the pre-thermalized stage (for initial conditions see [13 )• The decoupling temperature is set to Td cc ~ 130 MeV. In 
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the longitudinal invariance scenario, the temperature of the medium at any space-time point is defined by the proper 
time and the radial position, such that T — T(r, xj_). The flow velocity for any space-time rapidity value can be 
extracted from the flow at rj = 0. Indeed, for (3 = (/3 ± ,/3 z ), we have 

(3 ± (t, V ,x ± ) = /3±(T ' C ^°' XX) , ^(r,7 7 ,x ± )=tanhr ? . (26) 

The calculated photon spectra, and their different components, are shown in Fig. [4j The data are for Au-Au 
collisions, at the top RHIC energy, for two different centrality classes. The — 10% and — 20% classes are shown 
in the left and right panel, respectively. Note the data in the larger class extends to lower pt than that for the 
more central class, owing to a different experimental extraction technique [26]. Considering first that figure (left 
panel), the different contributions highlighted are those from hard primordial scatterings (prompt), which include 
the photons from Compton and annihilation events, together wit those from the fragmentation of jets. The photon 
spectrum associated with the interaction of jets with the thermal components of the quark-gluon plasma is labeled 
jet-QGP. The radiation from the thermal components of the quark-gluon plasma is shown, together with that from 
the thermal components of the hotgas of composite hadrons. The sum of the different contributions is the solid 
curve; the data are from PHENIX [2y]. It is seen that, for the physical case under consideration here, the jet-plasma 




FIG. 4: (Color online) Yield of photons in Au+Au collisions at RHIC, for two centrality classes: 0-20% (left panel) and 0-10% 
(right panel). The different elements of the theoretical calculation are described in the text. The data are from Refs. [2(|, and 
[29l | . respectively. 

photons are important to the theoretical interpretation of the experimental data in the window 2 < pt < 4 GeV. For 
smaller values of px, the emission from thermal media (whether QGP or hadron gas) represents a sizeable source. 
For the higher transverse momentum data, the radiation from hard collisions gradually take over the whole spectrum. 
This picture receives additional support from the higher p T data in the right panel. Most of that data is dominated 
by Compton and annihilation contributions calculated from pQCD. The jet-plasma sources are demanded only by 
the first two data points. The purely thermal contributions are subdominant in the entire range spanned by this 
figure. The fragmentation contribution to the real photon spectrum is small, owing mainly to the energy lost by the 
propagating jets. Here again, adding all of the sources produces a signal in agreement with the measured data. 

Another useful representation of the experimental data and a quantitative measure of the nuclear effects is provided 
by a plot of R AA (for real photons), shown in Fig. [5] The experimental data seem to show an interesting trend 
pointing towards diminishing values of R AA as pr grows. The experimental error bars are too large to permit a 
precise quantitative assessment, but different interesting possibilities and combinations may be considered. The 
left panel of Fig. [5] shows R AA calculated under different assumptions, and basically shows the importance of cold 
nuclear effects. The dashed curve shows the effect of the nuclear environment (shadowing) on the parton distribution 
function (pdf), while neglecting the specific isospin composition of the colliding nuclei. The full curve includes both 
isospin and shadowing contributions. The results of both calculations are systematically higher than the experimental 
data centroids, and exhibit a smaller slope than the one seen in the measurements, although the isospin effect can 
cause a 20% reduction at high-pr, as also found in Ref. [30| . The right panel includes medium effects, calculated 
as described earlier in the text; all curves except one contain jet-plasma photons, together with leading parton 
energy loss as evaluated with AMY. The dashed-dotted line shows the effect of neglecting the isospin content of 
the parton distribution functions. The double-dash dotted curve shows the scale-dependence of R AA , with the 
result of using Q — pt for the prompt contribution instead of Q = Pt/ V* used elsewhere in this work. The full 
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FIG. 5: (Color online) Nuclear modification factor of direct photon in central Au+Au collisions at RHIC in 2D+1 hydro, with 
a scale Q = pr/V2 in the prompt contribution. Left panel: effect of shadowing and isospin on the prompt contribution without 
medium effects. Righ panel: the effect of QGP and the scale is studied. The effect of a scale Q = pr is shown by the double 
dash-dotted line, while the effect of removing all photons produced from jet-medium interactions is shown by the dashed line. 
The result obtained without isospin effects is shown by the dot-dashed line. Data points are from PHENIX [29| . 



curve shows the nuclear modification factor evaluated with all sources described in this paper, together with the 
relativistic hydrodynamics evolution. Recall that the relativistic hydrodynamics modeling is constrained by a set of 
soft hadronic data [ill ]. The larger visible effect on the nuclear modification factor appears when jet-plasma photons 
are neglected (dashed line), causing a 30% reduction at pr = 8 GeV. The jets are however allowed to loose energy 
before fragmentation (like all cases in this panel). Because of the large errors, the data does not currently permit 
to choose between the cases where the jet-plasma photons are present or absent. However, it is important to realize 
that RZ A < 1 at higher values of pr, is a direct consequences of the fragmentation photons being affected by the 
energy loss of the fragmenting jet, as well as isospin effect in the nucleus-pdf. Should this trend, apparent in Figure 
[5J be confirmed experimentally, a quantitative link would exist between the high momentum nuclear modification 
factor of photons, and that of strongly interacting particles also born out of jet fragmentation. It is important for the 
same approach to reproduce both observables. Also, the large values of R\ A observed at pt < 6GeV/c (right panel 
of Fig. O are directly attributable to thermally-induced channels, in our approach. Our calculated results appear 
to overestimate the central values of the measured quantities (note however that the denominator of R AA is slightly 
underestimated at low pt by pQCD: correcting this will make our result correspondingly smaller), but smaller error 
bars would go a long in quantifying the medium-related processes. 
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FIG. 6: (Color online) Azimuthal anisotropy of direct photons in 20-40% central collisions at RHIC, within a 2D+1 hydro 
model. Dashed line : jet-plasma contributions; dot-dashed line: jet-fragmentation contribution; double dot-dashed line: thermal 
radiation of QGP; solid line: sum of QGP, prompt and hadronic gas contributions. The data are from Ref. (3lJ. 

We turn now to calculations and measurements of photon azimuthal anisotropy. This was discussed for low px 
photons in Ref. [28j], and for high pt photons in Ref. y|; both regions are treated here. Using Eq. flUJ) , (for real 
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photons) can be calculated for the different ingredients of the theoretical treatment, and compared to experimental 
data. This is done in Fig. [5] One observes that the net anisotropy (full line) is very low and in fact essentially 
zero for px > 4 GeV. For smaller values of the transverse momentum, our results are smaller than the central values 
of the experimental data (even if the error bars are large). The vi coefficient for the photons originating from the 
jet-plasma interactions is indeed negative [5[, but its magnitude is numerically smaller than in these earlier estimates. 
The explanation for this difference is two-fold. First, the present calculations relies on a realistic 2D+1 modeling [ll[ 
of the space-time evolution of the hot and dense medium, as opposed to using a simpler ID Bjorken expansion. The 
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FIG. 7: (Color online) The effect of flow and of the geometry of the initial jet-source profile, on the azimuthal anisotropy (left 
panel) and on the spectrum of photons (right panel). 

second reason has to do with the extreme sensitivity of on the initial conditions of the cooling and expanding source. 
This is seen in Figure[7] The geometry of the initial jet profile {Tab) is calculated here using a realistic Woods-Saxon 
distribution, unless noted otherwise. In the alternate case, thick spheres were used to calculate the nuclear overlap. 
A general observation is that is essentially flat as a function of pr, for all cases studied here. Considering first 
the jet-fragmentation contribution (in other words, the fragmentation contribution of the prompt production), the 
azimuthal anisotropy is positive and larger for the thick sphere geometry, and no radial flow. In the no-flow exercise, 
j3 is set to and the temperature at each point of the medium evolves from the initial profile Ti(r±) according to 
a Bjorken one-dimensional expansion T(rj_) = T i (rj^)(rj/r) 1 / 3 . The Woods-Saxon distribution, together with the 
flow bring down to a level of roughly 2%. The situation is similar for the jet-plasma component, the azimuthal 
anisotropy has a similar magnitude, but with an overall negative sign. The consequences of ignoring the radial flow 
and of varying the geometrical profile arc negligible for the spectrum of jet-plasma photons, as seen from the right 
panel of Figure [7J While the high-p^ QGP thermal photons, produced early in the medium evolution, are not much 
affected by the presence of transverse flow, they are largely affected (by a factor 3) by the choice of the geometrical 
profile which determine the temperature profile in the transverse plane, to which the thermal production rates are 
sensitive. At low-pr, thermal radiations dominate the spectrum and bring the total photon anisotropy coefficient up 
to 3-4% (FigO. This behaviour is however difficult to validate with the experimental data, owing again to the large 
size of the error bars. 

In summary, the azimuthal anisotropy is smaller (closer to zero) using a 2D+1 hydrodynamic model with realistic 
geometry than using a ID Bjorken expansion with a thick sphere initial jet distribution. The flow dynamics and the 
Woods-Saxon profile conspire to create a smaller geometrical anisotropy as seen by the traveling jets. However, these 
results do show that the photon v-i is rather sensitive to early time dynamics in relativistic nuclear collisions: precise 
measurements have the potential to stringently constrain evolution approaches. 

Lastly, similar calculations as the ones performed here for photons may be done for high momentum lepton pairs, 
using the techniques of Ref. [32^ . Even though no experimental measurements are yet available for the set of 
kinematical cuts we apply, it is nevertheless instructive to examine the relative importance and behavior of the different 
contributions shown in Figure [5J The components involving a plasma contribution are numerically important in the 
region px (GeV) < 8. These include the radiation from the thermal QCD plasma, the effect of jet plasma interactions 
(with and without the collinear enhancement germane to the many-body treatment). At higher transverse momenta, 
the spectrum shown here is taken over by the combination of prompt and fragmentation dileptons. We do not show 
here the pairs coming from the correlated semi-leptonic decay of heavy-quark mesons [32| . 
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FIG. 8: (Color online) Yield of dileptons in central 0-10% collisions at RHIC, within a 2D+1 hydro model. See main text for 
details. 



Assuming that the early phases of the relativistic collisions of heavy nuclei produce a quark-gluon plasma which 
later hadronizes into a hot gas of strongly interacting particles, we have calculated the electromagnetic signature 
of the different phases of the hot and dense dense, and evaluated the integrated signal. Within present error bars, 
we find agreement with the preliminary measurements of the PHENIX collaboration for the photon spectra, the 
photon nuclear modification factor, and the photon azimuthal anisotropy. Importantly, the dynamical evolution of 
the strongly interacting system is governed by relativistic hydrodynamics with no parameters in addition to the ones 
needed to quantitatively reproduce a large set of soft hadronic observables. Induced contributions from the QGP 
influence significantly the spectra, R AA , and elliptic flow of photons for px < 5 GeV. Our calculations predict larger 
values for R\ A and smaller values for v"l than suggested by the central values of the experimental data presently 
available [261 ] - The latter are, however still subject to considerable experimental uncertainties. At higher pt, our 
results confirm that electromagnetic observables are precise probes of early time dynamics, especially the signal's 
azimuthal anisotropy in momentum space. According to our calculations, the jet-plasma contributions dominate the 
photons yield in the window 2< pr <4 GeV. More precise value about the nuclear modification factor of photons 
would have the potentiel to discriminate between an inclusive and a non-inclusive jet-plasma contributions scenario. 
The extension of our findings to LHC energies is under way. 
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